Nonequilibrium Mesoscopic Conductors Driven by Reservoirs 
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In order to specify a nonequilibrium steady si^te of a quantum wire (QWR), one must connect 
reservoirs to it. Since reservoirs should be large 2d or 3d systems, the total system is a large and 
inhomogeneous 2d or 3d system, in which e-e interactions have the same strength in all regions. 
However, most theories of interacting electrons in QWR considered simplified Id models, in which 
reservoirs are absent or replaced with noninteracting Id leads. We first discuss fundamental problems 
of such theories in view of nonequilibrium statistical mechanics. We then present formulations which 
are free from such difficulties, and discuss what is going on in mesoscopic systems in nonequilibrium 
steady state. In particular, we point out important roles of energy corrections and non-mechanical 
forces, which are induced by a finite current. 



I. INTRODUCTION 

According to nonequilibrium thermodynamics, one can specify nonequilibrium states of macroscopic systems by 
specifying local values of thermodynamical quantities, such as the local density and the local temperature, because 
of the local equilibrium When one studies transport properties of a mesoscopic conductor (quantum wire 

(QWR)), however, the local equilibrium is not realized in it, because it is too small. Hence, in order to specify 
its nonequilibrium state uniquely, one must connect reservoirs to it, and specify their chemical potentials (/xl, /xr) 
instead of specifying the local quantities of the conductor (Fig. |l]). The reservoirs should be large (macroscopic) 2d 
or 3d systems. Therefore, to really understand transport properties, we must analyze such a composite system of the 
QWR and the 2d or 3d reservoirs. Although the QWR itself may be a homogeneous Id system, the total system is a 
2d or 3d inhomogeneous system without the translational symmetry. Moreover, many-body interactions are important 
both in the conductor and in the reservoirs: If electrons were free in a reservoir, neither could electrons be injected 
(absorbed) into (from) the conductor, nor could they relax to achieve the local equilibrium. However, most theories 
considered simplified Id models, in which reservoirs are absent or replaced with noninteracting Id leads |^-|l^. 
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FIG. 1. A two-terminal conductor composed of a QWR and reservoirs. 

In this paper, we study transport properties of a composite system of a QWR plus reservoirs, where e-e interactions 
are present in all regions. By critically reviewing theories of the conductance, we first point out fundamental problems 
of the theories in view of nonequilibrium statistical mechanics. We then present formulations which are free from 
such difficulties, and discuss what is going on in mesoscopic systems in nonequilibrium steady state. In particular, we 
point out important roles of energy corrections and non- mechanical forces, which are induced by a finite current. 
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II. A CRITICAL REVIEW OF THEORIES OF THE DC CONDUCTANCE 



In this section, we critically review theories of the DC conductance G of interacting electrons in a QWR. Note that 
two theories which predict different nonequilibrium states can (be adjusted to) give the same value of G (to agree with 
experiment). Hence, the comparison of the values of G among different theories is not sufficient. For definiteness, 
we consider a two-terminal conductor composed of a quantum wire (QWR) and two reservoirs (Fig. which are 
defined by a confining potential u'^, at zero temperature. Throughout this paper, we assume that u'^ is smooth and 
slowly-varying, so that electrons are not reflected by u'^ (i.e., the wavefunction evolves adiabatically) . We also assume 
that only the lowest subband of the QWR is occupied by electrons. A finite current / is induced by applying a finite 
difference Afi = fJ-h — /iR of chemical potentials between the two reservoirs, and the DC conductance is defined by 
G = {I)/{Afi/e) [|l^, where (/) is the average value of /. 

Let us consider a clean QWR, which has no impurities or defects. For non-interacting electrons the Landauer- 
Biittiker formula gives G = e'^/irh whereas G for interacting electrons has been a subject of controversy [p^ . 
Most theories before 1995 predicted that G should be "renormalized" by the e-e interactions as G = Kpe'^/nh, 
where Kp is a parameter characterizing the Tomonaga-Luttinger liquid (TLL) |p^|~|l9||. However, Tarucha et al. found 
experimentally that G ~ e^/nh for a QWR of Kp ~ 0. 7 [pOl . Then, several theoretical papers have been published to 
explain the absence of the renormalization of G P- |l^j2ll |~Although they concluded the same result, G = e'^/wh, the 
theoretical frameworks and the physics are very different from each other. Since most theories are based either on 
the Kubo formula (or, similar ones based on the adiabatic switching of an "external" field), or on the scattering 
theory, we review these two types of theories critically in this section. 

A. Problems and limitations of the Kubo formula when it is applied to mesoscopic conductors 

When one considers a physical system, it always interacts with other systems, Ri, R2, ■ • ■, which are called heat 
baths or reservoirs. Nonequilibrium properties of the system can be calculated if one knows the reduced density 
matrix ( = TrRi+R2+ - [Ctotai]- Here, Ctotai is the density operator of the total system, and TrFi,i+R2+ -- denotes the 
trace operation over reservoirs' degrees of freedom. To find C, Kubo p2[ assumed that the system is initially in its 
equilibrium state. Then an "external field" Eoxt is applied adiabatically(i.e., Eoxt oc e"*^'*'), which is a fictitious field 
because it does not always have its physical correspondence (see below). The time evolution of C was calculated using 
the von Neumann equation of an isolated system; i.e., it was assumed that the system were isolated from the reservoirs 
during the time evolution Because of these two assumptions (the fictitious field and isolated system), some 
conditions are required to get correct results by the Kubo formula. To examine the conditions, we must distinguish 
between non-dissipative responses (such as the DC magnetic susceptibility) and dissipative responses (such as the DC 
conductivity a). The non-dissipative responses are essentially equilibrium properties of the system; in fact, they can 
be calculated from equilibrium statistical mechanics. 

For non-dissipative responses, Kubo |^,^ and Suzuki |Q established the conditions for the validity of the Kubo 
formula, by comparing the formula with the results of equilibrium statistical mechanics: (i) The proper order should 
be taken in the limiting procedures of w, g — *■ and V —^ 00, where lo and q are the frequency and wavenumber of the 
external field, and V denotes the system volume, (ii) The dynamics of the system should have the following property; 

lim (iB(i))eq = (i)cq(S)eq, (1) 
t — >oo 

where (• • •)cq denotes the expectation value in the thermal equilibrium, and A and 13 are the operators whose correla- 
tion is evaluated in the Kubo formula. Any integrable models do not have this property [ P4|j2q -Pq]. Hence, the Kubo 
formula is not applicable to integrable models, such as the Luttinger model, even for (the simple case of) non-dissipative 
responses Q. 
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FIG. 2. Schematic plots of the chemical potential fi [M and the electrostatic potential 0, for (a) a macroscopic inhomogeneous 
conductor and (b) a mesoscopic conductor. For case (a), the local equilibrium is established, and thus fj, and (f) can be defined 
in all regions. The differences eA^ and A/i are not equal if one takes the differences between both ends of the conductor, 
whereas e A(/) = Afi if the differences are taken between the leads. For case (b) , /i cannot be defined in the QWR and boundary 
regions (although in some cases /j, could be defined separately for left- and right-going electrons), whereas can be defined in 
all regions. Similarly to case (a), eA(j) ^ Afi if one takes the differences between both ends of the QWR, whereas eA<j} — A/i if 
the differences are taken between the reservoirs. 



For dissipative responses, the conditions for the applicability of the Kubo formula would be stronger. Unfortunately, 
however, they are not completely clarified, and we here list some of known or suggested conditions for a: 
(i') Like as condition (i), the proper order should be taken in the limiting procedures. For a the order should be |2^] 

a= lim lim lim crformuia(g, w; 1^)- (2) 

(ii') Concerning condition (ii), a stronger condition seems necessary for dissipative responses: The closed system that 
is taken in the calculation of the Kubo formula should have the thermodynamical stability, i.e., it approaches the 
thermal equilibrium when it is initially subject to a macroscopic perturbation. (Otherwise, it would be unlikely for 
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the system to approach the correct steady state in the presence of an external field.) In classical Hamiltonian systems, 



this condition is almost equivalent to the "mixing property" 26-2^, which states that Eq. (|1|) should hold for any A 
and B, where (• • •)eq is now taken as the average over the equi-energy surface. It is this condition, rather than the 
"ergodicity" , that guarantees the thermodynamical stability p6|-^ . Although real physical systems should always 
have this property, some theoretical models do not. In particular, any integrable models do not have this property 

MM- 

(iii') We here suggest that all driving forces, including non- mechanical ones, should be identified |Q. In fact, the 
formula gives the current density in the following form, 

(J) = CformulaEext , (3) 

whereas the exact definition of a is given by noncquilibrium thermodynamics as jl]j 

(J) = -aV(/Ve) - L12V/3 = aE - aV(A*c/e) - ii2V/3. (4) 

Here, /3 denotes the inverse temperature, /i is the "chemical potential" which consists of a chemical portion fi^ and 
the electrostatic potential (j) p],pO|; 

fi = fic + ecf) (hence, A/i — A/ic + eA(j> for differences). (5) 

Hence, to evaluate a, one must find the relation between Eext and E, V fJ-c and V/3. In homogeneous systems, it is 
expected that Wfic — V/3 = 0, hence it is sufficient to find the relation between the fictitious field Ecxt and the real 
field E [| 10 3^. In inhomogeneous systems, however, V/Xc 7^ and/or V/3 7^ in gener al 13^ , as shown in Fig. ^ (a). 



Therefore, one must find the relation between Ecxt and these "non-mechanical forces" |29|]35|. (See section |y|.) 

Unfortunately, these conditions are not satisfied in theories based on simplified models of mesoscopic systems. 
For example, the Luttinger model 10 used in much literature does not satisfy conditions (i) and (ii) because it 
is integrable. To get reasonable results, subtle procedures, which have not been justified yet, were taken in actual 
calculations. Moreover, the non-mechanical forces have not been examined, although they would be important because 
a mesoscopic conductor (a QWR plus reservoirs) is an inhomogeneous system. 

We also mention limitation of the Kubo formula: it cannot be applied to the noncquilibrium noise (NEN) , which is 
the current fluctuation in the presence of a finite current (/) (= GAfi) The NEN at low frequency, ((5/^)'^-°, 

is usually proportional to |(/)| cx |A/i|. However, the Kubo formula assumes power series expansion about A/i ~ 0, 
hence cannot give any function of |A/i| p9|,^. 



In sections [II and [V, we present other formulations which are free from these problems and limitations. These 
formulations clarify what is going on in noncquilibrium mesoscopic conductors, because one can find the noncquilibrium 
steady state. This is impossible by the Kubo formula because it evaluates correlation functions in the equilibrium 
state. 



B. Scattering-theoretical approaches 

In view of many problems and limitations of the Kubo formula, it is natural to try to generalize Landauer's theory 
to treat conductors with many-body interactions. Namely, the DC conductance may be given in terms of the 
scattering matrix (S matrix) for interacting electrons ||7|, pT| , ^ . 

The advantages of the scattering-theoretical approaches may be as follows: (i) Neither the translation of A^oxt into 
Afi nor the subtle limiting procedures of w, q and V is necessary, (ii) There is no need for the mixing property of the 
Id Hamiltonian Hi. Hence, Hi can be the Hamiltonian of integrable Id systems such as the TLL. (iii) In contrast to 
the Kubo formula, one can calculate the NEN |@,|^-|6) . 

However, to define the S matrix, one must define incoming and outgoing states. Although they can be defined 
trivially for free electrons, it is nontrivial in the presence of many-body interactions. In high-energy physics, they 
are defined based on the asymptotic condition, which assumes that particles behave like free (but renormalized) ones 
as i ^ ±00, i.e., before and after the collision |Q. For example, an electron (in the vacuum) before or after the 
collision becomes a localized "cloud" of electrons and positrons, which extend only over the Compton length, and this 
cloud can be regarded as a renormalized electron. In condensed-matter physics, on the other hand, the asymptotic 
condition is not satisfied for electrons in metals and doped semiconductors. In fact, elementary excitations (Landau's 
quasi particles) are accompanied with the backflow, which extends all over the crystal in contradiction to the 
asymptotic condition. Because of this fundamental difficulty, the scattering approaches to mesoscopic conductors 
replaced the reservoirs with Id leads in which electrons are free Therefore, real reservoirs, in which electrons 

behave as 2d or 3d interacting electrons, have not been treated by the scattering-theoretical approaches. 
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III. COMBINED USE OF MICROSCOPIC THEORY AND THERMODYNAMICS 



The basic idea of this method ||2l| is as foUows: Since a QWR is a small system, and is most important, it 
should be treated with a full quantum theory. On the other hand, reservoirs are large systems whose dynamics is 
complicated, hence it could be treated with thermodynamics (in a wide sense). Utilizing these observations, we shall 
develop thermodynamical arguments to find the nonequilibrium steady state that is realized when a finite is 
applied between the reservoirs. This is the key of this method because when the steady state is found, G (and other 
observables) can be calculated by straightforward calculations. Although in some cases formal calculations can be 
performed without finding the steady state we stress that such formal theories are incomplete because another 
theory is required to relate A/i of such theories with A/i of the reservoirs, by which G is defined. 

An advantage of the present method is that we do not need to find the relation between A(/)oxt and A/j, because (/) 
is directly calculated as a function of A//. Another advantage is that it is applicable to NEN and nonlinear responses 
because nonequilibrium steady state is directly obtained. 



A. Conductance of the Id Fermi liquid 



It is generally believed that a Id interacting electron system is not the Fermi liquid (F L) but the Tomonaga- 
Luttinger hquid (TLL) p6|-p^. For this reason, many papers on Id systems [|3|-|6|,p|,p Hll]jl2| ] use the word FL to 
indicate non-interacting electrons, i.e., a Fermi gas. However, we do not use such a misleading terminology; by a FL 
we mean interacting quasi-particles. Since the backflow is induced by the interaction pH , the Landauer's argument 
of non-interacting particles cannot be applied to a FL. On the other hand, real systems have finite length and 
finite intersubband energies, in contradiction to the assumptions of the TLL. Hence, some real systems might be well 
described as a FL. Therefore, G of a FL is non-trivial and interesting jl^]. Furthermore, we will show in section 
^ that the results for the FL suggest very important phenomena that is characteristic to nonequilibrium states 
of inhomogeneous systems. Note also that the following calculations look similar to the derivation of fundamental 
relations in the theory of the FL |4l]. However, G of mesoscopic conductors was not calculated in such calculations. 
The most important point to evaluate G is to find the nonequilibrium steady state. 

We find the nonequilibrium steady state using a thermodynamical argument as follows: In the reservoirs, electrons 
behave as a 2d or 3d (depending on the thickness of the reservoir regions) FL. Since we have assumed that u'^ is smooth 
and slowly-varying, a 2d or 3d quasi-particle in a reservoir, together with its backflow, can evolve adiabatically into a 
Id quasi-particle and its backflow in the QWR, without reflection. In this adiabatic evolution, the quasi-particle mass 
m* and the Landau parameters / also evolve adiabatically, and the energy is conserved. Therefore, quasi-particles 
with e(fc > 0) < Ml are injected from the left reservoir. Here, e is the quasi-particle energy; 



e{k) 



2m* 



(6) 



where Sn{k) = n{k) — 6(|fc| < kp), with n{k) being the quasi-particle distribution. The last term of this expression 
represents energy correction by interactions among quasi-particles p| . On the other hand, a quasi-hole below 
should not be injected because otherwise the recombination of a quasi-particle with the quasi-hole would produce 
excess entropy, in contradiction with the principle of minimum entropy production. Similarly, quasi-particles with 
e{k < 0) < /iR are injected from the right reservoir, with no quasi-holes are injected below /ir. Therefore, the 
nonequilibrium steady state under a finite A/i = fJ-L — ^J'R should be the "shifted Fermi state" , in which quasi-particle 
states with e{k > 0) < /il and e{k < 0) < fiR are all occupied. Hence, the right- (left-) going quasi-particles have the 
chemical potential = /iL (m+ = fJ-n)- Considering also the charge neutrality, we can write the distribution function 
as 



Then, Eq. (|) yields 



where /++ = /(fcp, kp) and /+_ 



n{k) = e(|fc- g| < fcp). 



n'kl 

2m* 



± hq 



hkp 
m* 



/++ ~ /-+ 
2tt 



f{kp,—kp)~f{—kp,kp). Hence, 

' ^kp ^ /++ - /+- 
TO* 27r 



A/i = 2fiq 



(7) 



(8) 



(9) 
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On the other hand, considering the spin degeneracy, (/) is calculated as 



(/) = 2e 



/++ ~ /-+ 
27r 



Here, the /+±-dependent terms represent the backflow. Since the same factor appears in Eq. 
conductance is independent of m* and /+±; 



G 



(!) 



e 



(10) 

wc find that the 
(11) 



Since we have identified the nonequilibrium steady state, we can calculate not only G but also other nonequilibrium 
properties such as the NEN 

It is instructive to represent Eqs. in terms of the bare parameters. As in the case of 3d Fermi liquid 

we can show that Ml 



Tik-p 
m 



Tik-p 
m* 



/++ - h 
2tt 



(12) 



Hence, we can rewrite Eq. (10) as (/) = 2e{q/T:){hk-p /m). Therefore, quasi particles (whose group velocity is hkp/m*) 
plus their backflows carry exactly the same current as the bare particles, for the same q, i.e., for the same shifted 
Fermi distribution. On the other hand, Eq. (||) is rewritten as fi± = h^k^/2m* ± h^qk-p/m. Although fi± ^ of 
bare particles], A/x — [A/it of bare particles] for the same q. These facts result in the independence of G on the Landau 
parameters. 



B. Conductance of the Tomonaga-Luttinger liquid 

We now consider a clean TLL . The low-energy dynamics of a TLL is described by the charge (p) and spin (cr) 
excitations (whose quantum numbers are and iV^, respectively, where q ^ Q denotes the wavenumber), and the 
zero modes (quantum numbers A^^, A^±) Il6 -19|. The eigenenergy is given by 

^ = E «''E^i'?K + iJ ^ bXr(A^;+A^-)'+«j(A^;-Ar-)'], (13) 

v—p^a q u—p^a 

where v'^ = v'^/K^ and = K^v'^ [v = p,a). Here, the parameters v'^ and are renormalized by the e-e 
interactions (except that = 1 by the SU(2) symmetry). The DC current is given by 

(I) ^2ev>;j{N^~ N^)/C. (14) 

We apply a thermodynamical argument to find the nonequilibrium steady state. Unlike the FL case, there is no 
adiabatic continuity between the TLL in the QWR and the FL in the reservoirs. We therefore argue differently: In 
the linear response regime the steady state must be the state with the minimum energy among states which satisfy 
given external conditions. Otherwise, the system would be unstable and would evolve into a state with lower energy. 
For our purpose, it is convenient to take the value of (/) as the given external condition. Then, from Eqs. (^3|) and (|lj), 
we find that the steady state should be the state with NP = = (for all q), A^^ + A;^ = 0, A^^+A^^ = N^-NZ = 0, 
and A'^^ — A'^^ > 0. This state may be called the "shifted Fermi state" of the TLL. Furthermore, in the steady state, 
electrons in the left reservoir and right-going electrons in the TLL should be in the "chemical equilibrium" , in which 
electrons in the FL phase are transformed into right-going electrons in the TLL phase at a constant rate. Therefore, 
their chemical potentials should be equal 

= M^t,- ^ - + ± v^jiN!l - A^)], (15) 



where we have used Eq. (|l^). Hence, 



AM = ML-m = ^^^5W-^-)- (16) 
By dividing Eq. ( p^ by this expression, we obtain the same result for G as Eq. (pi]), in agreement with experiment 



Since we have identified the nonequilibrium steady state, we can calculate not only G but also other nonequilibrium 
properties such as the NEN plf . 
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IV. PROJECTION THEORY 



Although we have successfully found the nonequilibrium steady state of interacting electrons in section III, a 
possible objection against the formulation may be that the theory is rather intuitive. In this section, we present 
a full statistical-mechanical theory Q,^, which is free from such an objection. In this theory, we start from the 
Hamiltonian of 3d interacting electrons confined in the composite system of the QWR and reservoirs. Fig. I. This 
original system is projected onto an effective Id system, and the equation of motion for the reduced density operator 
of the Id system is derived. From this equation, we can find the nonequilibrium steady state as a function of A/x 
between the reservoirs. This allows us to evaluate various nonequilibrium properties. 



A. Decomposition of the 3d electron field 



We start from the 3d electron field ■(/;(r) subject to a confining potential u'^{r) (which defines the QWR and two 
reservoirs connected to it), impurity potential M'(r) (whose average u' is absorbed in u'^(r), hence u' = 0), external 
electrostatic potential 0cxt(r), and the e-e interaction of equal strength v{r — r') in all regions |Q: 

^2 



— + u^iv) + u\t) + e0ext(r) 
Im 



V'(r) 



-I — / (Pr / (Fr' p{y)v{y -y')p{y'), 



(17) 



where p(r) is the charge density. We will find the nonequilibrium steady state for A/i > 0. For this state, (p(r)) ^ 0, 
which gives rise to a long-range force. We extract it as the rcnormalization of the electrostatic potential 



e0(r) = e0cxt(r) + / v{y - r')(p(r')) 



(18) 



and a c-number Kiv Namely, H is recast in terms of 5p{y) = p{y) — {p{y)) as 



H 



(fr V3^(r) 



-— + M'=(r) + u\y) + e0(r) 
Zra 



V'(r) 



+ - / (fr / (Pr' 5p{y)v{y - y')6p{y') + Va, 



(19) 



To decompose V'(r), we consider the single-body part of H . Recall that u'^{y) is assumed to be smooth and slowly- 
varying, to avoid undesirable reflections at the QWR-reservoir boundaries. In this case, the single-body Schrodinger 
equation 

\{Y)^e^[Y) (20) 



+ u''{y) + u\y) + e0(r) 

2m 



has solutions that propagate through the QWR having the energy e ~ e-p |4^] ; 



Kk{x)dx 



(21) 



Here, ip'^{y, z;x) is the wavefunction of the lowest subband at x, representing the confinement in the lateral (yz) 
directions, and C is the normalization length in the x direction. All the other modes are denoted by </?iy(r), which 
includes solutions that are localized in either reservoir, and extended solutions whose e are not close to ep- Since any 
fmiction of r can be expanded in terms of (^fc(r)'s and (py(r)'s, so is the r dependence of the electron field operator; 



= ^Ck(pk{r) + ^di,(pi,{Y) = f^{y,z;x)'ipi{x) +'4'r{y). 



(22) 



The 3d electron field has thus been decomposed into the Id field and the 3d field 'i/'R(r), which we call the "reservoir 
field." It can be decomposed into the low-energy components ip^,^ and tpn^ (which are localized in the left and right 
reservoirs, Rl and Rr, respectively) and the high-energy component ^-r.h 



■IpR = V'Rl + "V^Rr + "^Rh- 

For low-energy phenomena, we can take ■i/'R = "^Rl + V'Rr- 



(23) 
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B. Hamiltonian for the Id and the reservoir fields 



By expressing H in terms of ipi and tAr, we obtain the single-body part of ipi (denoted by Hi), the interaction 
(Vii), the V^i-V^R interaction (V^ir), the ipR-i^n interaction (Vrr), and the single-body part of i/jr (-^r) [Q- By the 
screening effect of V'rr, V^r is renormalized as the screened interaction Vir. Similarly, by the screening effect of Vir, 
is renormalized as the screened interaction Vn. We therefore recast V^i + t^R + ^rr Vn + Vir -I- Vrr, where 
ViR and Vrr have no screening effects on Vn and Vir . In this way, H is decomposed as 

H^Hi+ Fir + Hr, (24) 

where Hi = H^ + Vn is the Hamiltonian for "01(2;), ^^r = ^^r + Vrr is the one for V'R(r), and Vir is the interaction 
between ipi{x) and ip-[i{r). In particular. Hi is evaluated as 



Hi= j dx i:\{x) 



^^9^ ±r ^ ir ^ 
2m ox'^ 



■^i{x) 



dx I dx' Spi{x,t)vii(x,x')5pi(x' ,t). (25) 



1 

+ 2 

Here, is the subband energy [Q, Spi(x,t) = il'l{x)il'i{x) — (?/'|(a;)'0i(x)) is the density fluctuation of the Id field, 
and 

u[ix) = J J dydz |^^(y,z;a:)|2 u\r), (26) 
vii{x,x') = J J J Jdydzdy'dz' |^^(2/, z; x)| V^(r, r') \^^{y' , z';x')\^ (27) 

arc the impurity and two-body potentials for the Id field, where v^'^ is the screened two-body potential for tp. 

It is seen that u\{x) is the average of the random potential ^^(r) over the lateral wavefunction (p-^, which, as a 
function of y, is localized in a region of width ~ W{x) for each x. Here, W{x) denotes the width of the region in 
which electrons are confined (Fig. |l|). From these observations, we can show that 

u[{x) - u\r) for x e QWR, (28) 
|w']^(x)| oc l/^/W{x) for X e a reservoir. (29) 

In a similar manner, the two-body potential for the Id field behaves as 

vii{x,x') <v'"=i[\x-x'\^ + W{Of]^^^) for X or x' e QWR, (30) 
Vii{x,x') ~ {r^'^ /W)v^'^ , for x ^ x' £ a reservoir. (31) 

where r^'^ denotes the range of v^'^, and v^'^ the average of v^^ in the region |r — r'| < r*''^. We now assume that the 
width of the reservoirs is very large; 

W{x) ^ oo as X ^ ±oo. (32) 

Then, Eqs. (||)-(|l]) yield u\{x) as X ^ ±cxD, and wii(x,x') ^ as x or x' ^ ±oo. Namely, Hi represents 
interacting ^pi{x) field that gets free as x ^ ±oo. On the other hand, the interaction Vir between ipi{x) and ipR{r) 
becomes stronger in the reservoir regions, whereas it is negligible in the QWR because at low energies tpnir) does not 
penetrate into the QWR. Therefore, the Id field '0i(x) is subject to different scatterings in different regions of x: In 
the QWR, ipi is scattered by the ipi-ipi interaction and the impurity potentials, whereas ipi is excited and attenuated 
by the reservoir field in the reservoir regions through the ipi'^'R interactions (Fig. S). 
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FIG. 3. Schematic diagram of the strengths of scatterings of the Id field. 



C. Equation of motion for the reduced density operator 



We have successfully rewritten the Hamiltonian H in terms of t/ji and ip^ . We must go one step further because 
H = Hi + ViR + -ffR describes very complicated dynamics, and thus the von Neumann equation for the density 
operator 



Hi + ViR + Hn,at) 



(33) 



is impossible to solve. This unsolvability guarantees the thermodynamical stability ( mixing property) of the total system 



|p6|-p8| . We turn this fact to our own advantage, and reduce the theory to a tractable one 45|. The basic idea is 
as follows: Hi describes the Id correlated electrons, and is most important. Hence, it should be given a full quantum- 
mechanical treatment. Concerning T4r, on the other hand, multiple interactions by Vir seem unimportant. Hence, 
we may treat it by a second-order perturbation theory [Q. For iJ^, it describes the 2d or 3d interacting electrons, for 
which many properties are well known, and we can utilize the established results. Moreover, since the reservoirs are 
large, we can assume the local equilibrium: both reservoirs are in their equilibrium states with the chemical potentials 
/iL and ^R, respectively. We denote the reduced density operator of the reservoir field for this local equilibrium state 
by Cr- 

From these observations, we may project out the reservoir field ?Ar as follows. Consider the reduced density operator 
for the Id field; C,i{t) = TrR[^(i)]. Up to the second order in Vir Q, the equation of motion of in the interaction 
picture of Hi + 7?r , is evaluated as 



5 ; . ^ -1 



dt' 



'TrR([ViR(0,[PiR(i'),CRCiW 



(34) 



where we have used the fact that Ci(0 the interaction picture varies only slowly, so that Ci(^') — Ci(0 the 
correlation time of (ViR(i)ViR(i')) [Hi- This equation represents that C,i is driven by two reservoirs, which have 
different chemical potentials /iL and /iR, through the V'fV'R interaction. Since the trace is taken over the reservoir 
field, Eq. (Q) is a closed equation for and Its steady solution represents the nonequilibrium steady state of the 
Id field driven by the reservoirs. 



D. Current of the Id field 

We now turn to observables. We are most interested in the ioia/ current / which is given by /(x, t) = J J dydzJ^i^jt), 

where Jx denotes the x component of the current density. Note that / is different from the current of the Id field 
defined by || 
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/i(x,t) 



e 
2m 



Tpl{x,t)\-—^l{x,t) \+h.C. 



(35) 



The expectation values of J, /i, and /r (the current carried by the reservoir field) are schematically plotted in Fig. 
^, which shows that (J) is mainly carried by (/i) in the QWR and by (/r) in the reservoirs, respectively. The 
transformation between Ii and /r is caused by Vir, and thus Ii is not conserved: + ^Pi 7^ 0. At first sight, 

these facts might seem to cause difficulties in calculating (/) and {51"^) from C,i. Fortunately, however, we can show 
that for any (nonequilibrium) steady state (/) and ((5/^)'^-° are independent of a;, and that 



{51^ 



(I) ^ (h) at x ~ 
= (5/2)'^-0 at x ~ 



(36) 
(37) 



where a; = corresponds to the center of the QWR. Therefore, to calculate (/) and ((5/^), it is sufficient to calculate 
(/i) and ((5 7^)"-" at a; ~ 0, which can be calculated from Ci- Therefore, we have successfully reduced the 3d problem, 
Eq. (|l|), into the efllective Id problem, Eqs. (|2|), (g) and (|3|). 



<I> 



<Ir> 





Reservoir 



Quantum Wire 



Reservoir 



X 



FIG. 4. Schematic plots of the expectation values (J), (Ii), and (/r), of the currents carried by tp, ipi, and ^p-R., respectively. 

Actual calculations can be conveniently performed as follows. Although (/) and ((5/^) have both low- and high- 
frequency components, we are only interested in the low- frequency components, which are denoted by I{t) and Sl'^{t), 
respectively. They are given by 



m 



dt'f{t'-tmi{t')at% 



(38) 



and similarly for 5P{t). Here, fit' — t) is a filter function that is finite only in the region t — t/2 ^ t' < t + t/2, 
where l/r ~ the highest frequency of interest. From Eqs. (p^), (|3^)-(38), we can construct the equations for I{t) and 
51^{t). They can be solved more easily than the equation for Ci, and the solutions fully describe the low- frequency 
behaviors of (J) and {5P). 

In the following, we present the results for lit) for the case of impurity scatterings and for the the case of e-e 
interactions. 



E. Application of the projection theory to the case where impurity scatterings are present in all regions 



When electrons are scattered by impurities (one-body potentials) in all regions including reservoirs, whereas many- 
body scatterings are negligible, the '0i"'0R interaction is given by 



IR 



d^r'ii}\{x)Lp^*{y, z; x)u\y)iP^{y) + h.c. 



(39) 



If we put Ya{v) = (y5^*(2/, z; a;)M'(r)ViR^ (r) [a — L, R), then Eq. (p3) becomes 
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dt 



/ dt' dPr d\'\ (i;t(r,i)f„(r',t')) ^(x, i), (x', OCi W 



+ (r„(r,t)yi(r',t')) v}(a^,t),V'i(x',i')CiW 



h.c. 



(40) 



where (• • ■)a denotes the expectation value for the cquihbrium state of reservoir R^, which has the chemical potential 
^a- After careful calculations using Eqs. (^-(^, and considering the spin degeneracy, we find 



dt 



lit) = -7 [I{t) - /.tcady] , 



(41) 



where /steady = (e/7r?i)TA/i, and 7 ~ (27r/?i)nimp|w'pI'F- Here, n-mip is the impurity density, u' denotes the potential 
of an impurity (w'(r) = u'5(r — r^)), and Dp is the density of states per unit volume, I?(/^l) — T^if^n) = 2?f- It is 
seen that / approaches /steady as i ^ 00. Therefore, the DC conductance is given by 



G = T 



(42) 



in agreement with the Landauer-Biittiker formula | p^ . Moreover, we find that the steady state is stable: For any 
(small) deviation from the steady state, / relaxes to the value /steady, with the relaxation constant 7. 



F. Application of the projection theory to the case where e-e scatterings are present in all regions 

When the e-e interaction is important in all regions including reservoirs, whereas impurity scatterings are negligible, 
we find that the most relevant term of the 'ipi-4'K interaction is given by 

ViR = jd\ y"dVV'I(x)(^^*(2/,z;a;)V;R(r)«^^(r,r')V'^(r')^R(r')+h.c. (43) 

By this interaction, an electron is scattered into the QWR through the collision of two electrons in a reservoir, or, an 
electron in the QWR is absorbed in a reservoir. By putting Fa(r) = Jd^r'ip^*{y,z]x)'ipYi^{Y)v^'^{r,r')'il)\^^{Y')'ij>Yi^{Y') 

(for a = L, R), which differs from Ya of the impurity-scattering case, we obtain the equation of motion for ^1 in the 
same form as Eq. (^0|). To derive the equation of motion for / from that equation, we need to calculate the correlation 
functions in the reservoirs, {Ya{t)Y^[t'))a and (yj(t)r„(i'))a- We can easily calculate them using well-known results 
for the 2D or 3D FL because the reservoir electrons are believed to be the 2D or 3D FL. We also need correlation 
functions of the Id field. They are quite different depending on the nature (FL or TLL) of the electrons in the QWR. 

In the case where ^"1 behaves as a Id FL, we obtain the equation for I{t) in the same form as Eq. (^Tj), but now 7 
is a function of the e-e interaction parameters, and 

/steady = 4 E " [©(^L " e{k)) - e(MR ^ e{-k))]. (44) 

fc>o 

Here, m is the hare mass, Q is the step function, and e(±fc) denotes the Id quasi-particle energy, Eq. (^), in the shifted 
Fermi state. Note that if we simply took e(±fc) = ?i^fc^/2TO*, then /steady = [m* /m){e/Trh)Afj,, hence the conductance 
would be renormalized by the factor m*/m. However, the correct expression (^ shows that e(±A;) are modified in 
the presence of a finite current, and the correction terms are proportional to g oc /. As a result, the injection of an 
electron becomes easier or harder as compared with the case of / = 0. This automatically "calibrates" the number of 
injected electrons, and we obtain 



-'steady — ^ 

m 



m Ztt 



-1 



where we have used Eq. (|T^). Therefore, G = e'^/nh. Here, the interaction parameters of the Id field are canceled in 
G, and those of the reservoir field are absorbed in 7. These observations confirm the results of section^ the shifted 
Fermi state is realized as the noncquilibrium steady state, and the conductance is quantized. 

The application of the projection theory to the case where ipi behaves as a TLL will be a subject of future study. 
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G. Advantages of the projection theory 

A disadvantage of the projection theory is that calculations of G become rather hard as compared with the simple 
theories that are reviewed in section However, the simple theories have many problems and limitations, as discussed 
there. The projection theory is free from such problems and limitations, and has the following advantages: (i) The 
value of (/) for the nonequilibrium state is directly calculated as a function of A/j,. Hence, neither the translation 
of A0cxt into nor the subtle limiting procedures of w, q and V is necessary, (ii) There is no need for the mixing 
property of the Id Hamiltonian Hi. Hence, Hi can be the Hamiltonian of integrable Id systems such as the TLL. (iii) 
In contrast to the Kubo formula, which evaluates transport coefficients from equilibrium fluctuations, the projection 
theory gives the nonequilibrium steady state. This allows us to discuss what Id state is realized and how the current 
is injected from the reservoirs. Moreover, we can calculate the NEN and nonlinear responses, (iv) The projection 
theory can describe the relaxation to the nonequilibrium steady state. This allows us to study the stability and the 
relaxation time of the nonequilibrium state. 



V. APPEARANCE OF A NON-MECHANICAL FORCE 



We here discuss the applicability of the Kubo formula to inhomogeneous systems. The general conclusion of this 
section is independent of natures (such as a FL or TLL) and the dimensionality of the electron system. Hence, we 



will use the results for the Id FL, which are obtained in section III A and confirmed by a full statistical-mechanical 



theory in section IV F 



The original form of the Kubo formula gives a conductivity that corresponds to the following conductance; 
(/)/A0oxt = GKubo Izuyama suggested that the conductance should be {I)/A(j> = Gizuyama, by consider- 

ing the screening of 0cxt POl ; ^ • On the other hand, the exact definition of the conductance is G = (/) / A/i [|l|,|j . 
For macroscopic inhomogeneous conductors, eA(f) ^ in general if one takes the differences between both ends of 
the conductor, as sketched in Fig. §(a). Therefore, G ^ GKubo, Gizuyama in such a case. Hence, to obtain the correct 
value of G by the Kubo formula, one must find the relation between A0oxt and A/x. Unfortunately, no systematic 
way of doing this has been developed. 

The same can be said for mesoscopic conductors. Fig. |^(b), for which eAcf) ^ Afi in general if one takes the 
differences between both ends of the QWR. Therefore, G ^ GKubo, Gizuyama. It is only for fortunate cases that GKubo 
or Gizuyama coincidcs with G. For example, Kawabata calculated Gizuyama for the case where the backward 
scattering with amplitude V{2kF) is present, which had been neglected in the previous calculations. He found that 



r — 

^Izuyama 



^ , V{2kF) 



2nhvF 



(46) 



However, this result disagrees with G obtained in the previous sections. The origin of this discrepancy may be 
understood as follows. By taking the Fourier transforms of both sides of Eq. (|l^), we can see that only the g ~ 
component of the two-body potential v contributes to the screening of the electrostatic potential. On the other 

hand, both g ~ (forward) and q ~ 2fcF (backward) components of v contribute the Landau parameter , i.e., 

/+_ = /f_°™ard _,_ ybackward_ Therefore, Eq. (|) shows that A^ has a term (proportional to /backward-, ^^^^^ cannot be 
interpreted as coming from the screening of (j)cxt- If we interpret this term in terms of nonequilibrium thermodynamics 
(although it is not fully applicable because the local equilibrium is not established), the term may be interpreted as 
a non-mechanical force in Eq. (j|): 

A^e ^ -(?ig/7r)/f_^'^™-'i + • • • . (47) 

Here, • • • accounts for possible contributions from and/or jfo™ard gjj^f.g ^ (/), so is A/ic. This means that 
a finite current (/) induces a finite non- mechanical force A/ic, and (/) is driven by both eAcj) and A/ic in the steady 
state. Hence, G = (/)/(A/i/e) is not equal to either GKubo = (-^)/A0cxt or Gizuyama = {I) I ^4>- Therefore, the Kubo 
formula cannot give the correct value of G if the q ~ 2k-p component of the two-body potential is non-negligible, even 
if the screening of A0oxt is correctly taken into account, because a non-mechanical force is inevitably induced. Note 
that this is not the unique problem of the Kubo formula, but a common problem of many microscopic theories which 
calculate a nonequilibrium state by applying a mechanical force. 

A possible way of getting the correct result by the Kubo formula would be to apply the formula to a larger system 
that includes the homogeneous reservoirs or leads ||5^: in that case, eA0 = A/i, as sketched in Fig. ||, and thus 
G — Gizuyama- Howcvcr, this seems very difficult because it is almost equivalent to trying to solve the Schrodinger 
equation of the total system, including complicated processes that lead to the mixing property and to the equality 
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eA(^ = A/Lt. Another possible solution may be to apply Zubarev's method 0, which, to the authors' knowledge, 
has not been applied to interacting electrons in mesoscopic conductors. However, one must also include (a part 
of) reservoirs into the Hamiltonian because Zubarev's method assumes that macroscopic variables (such as /z) are 
well-defined in the nonequilibrium steady state. As compared with these approaches, the formulations presented in 
and would be simpler ways of getting correct results which include effects of non-mechanical forces. 



sections ID 



VI. DEVIATION FROM THE QUANTIZED CONDUCTANCE 



For a clean QWR, we have obtained the quantized value G = /ttTi in both cases of the FL and the TLL in 



sections [II A, IIIB and IV F, using different formulations. The essential assumptions leading to this result are the 
following, (i) The QWR is clean enough and the temperature is low enough (zero temperature has been assumed 
for simplicity), so that scatterings by impurities, defects or phonons are negligible and e-e interactions are the only 
scattering mechanism, (ii) The boundaries between the QWR and reservoirs are smooth and slowly-varying so that 
reflections at the boundaries are absent, (iii) The reservoirs are large enough, so that they remain at equilibrium even 
in the presence of a finite current between the reservoirs through the QWR. 

When some of these assumptions are not satisfied the observed conductance may deviate from the quantized 
value. For example, if boundary reflections are non-negligible, the transmittance T (calculated from the single-body 
Schrodinger equation) between the reservoirs through the QWR is reduced. This results in the reduction of G by 
the factor T for non-interacting electrons. For interacting electrons, G will be further reduced for the TLL because 
the TLL will be "pinned" by the reflection potential at the boundaries. This can be understood simply as follows: 
Although the TLL of infinite length is a liquid, for which a long-range order is absent, it behaves like a solid at a short 
distance. Hence, the TLL is pinned by a local potential, like a charge-density wave is. This is the physical origin 
of the vanishing G (at zero temperature) for the case where a potential barrier is located in the TLL UJ^. Since 
the pinning occurs irrespective of the position of the local potential, the TTL would be pinned also by the boundary 
reflecti ons. N ote that if one neglects the weakening of vn in reservoirs (due to the broadening of W{x), as shown in 
section IV B), the TTL would then be pinned also by impurities in reservoirs pO| . 

Another example is dissipation by, say, phonon emission. By the dissipation, the Id system will lose any correlations 
over a distance Lrix, where L^\^ is the "maximal energy relaxation length" which is generally longer than the 

simple dephasing length (over which an energy correlation may be able to survive). In such a case the Id system of 
length L (> Lrix) will behave as a series of independent conductors of length irix. One will then observe Ohm's law 



Gobs - (irlx/i) X (eV7r?i). 



(48) 
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